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Abstract. We describe the design of a Kalman filter that identifies suspension violin 
modes in an interferometric gravitational wave detectors data channel. We demonstrate 
the filter's effectiveness by applying it to data taken on the LIGO 40M prototype. 

I MOTIVATION 

The wire suspensions of the interferometer test masses are a conduit for environ- 
mental noise to enter the gravitational wave data channel; additionally, they are a 
source of noise themselves. Both the pendulum mode (whose frequency is out of 
the interferometers band) and the suspension wire violin modes (whose fundamen- 
tal mode frequencies are in-band) are energized by their contact with the thermal 
bath. Owing to their weak damping this energy is strongly concentrated about the 
mode resonant frequencies. At the fundamental violin mode frequencies the ther- 
mal noise dominates the other noise sources by approximately 50 dB. In addition 
to this thermal noise component, non-thermal excitations of the suspension wires 
(e.g., sudden creep events) can lead to excitations in the interferometer output. 

These narrow band features are instrumental artifacts, not gravitational waves: 
along with other artifacts they should be removed from the data before it is studied 
for the presence of other signals. As instrumental artifacts, however, they carry 
important information about the instrument's state. Thermal and technical noise 
that disturbs the suspension excites these modes and move the mirrors, leading to 
an artifact in the gravitational wave channel. Gravitational waves give rise to a 
signal by changing the distance between the mirrors, but do not move or otherwise 
excite the suspension modes. Correspondingly, if we can determine the mode state 
as a function of time we have a way of eliminating a broad class of technical noise 
sources that might otherwise masquerade as gravitational wave bursts. 

How can we identify the state of the suspension violin modes, given only the 
gravitational wave channel? This is a classic problem in data analysis, generally 
addressed by a Kalman filter. We have developed such a filter for use in LIGO and 



shown, using data taken in November 1994 at the LIGO 40M prototype, that it 
enables us to follow the state of the violin modes independently of the other noise 
sources — both technical and fundamental — that contribute to the gravitational 
wave channel. 



II THE KALMAN FILTER 



The Kalman Filter [1] is a mechanism for predicting the multi-dimensional state 
of a dynamical system from a multi-dimensional observable. The system is assumed 
to evolve linearly and the observable is assumed to be linearly related to the state, 
denoting the system state x we have (for discrete time series): 

x[fc]=A-x[ib-l]+w[A;-l]. (1) 

We assume that the system is driven by a stochastic force, referred to as the process 
noise and denoted w. The state dynamics determine the linear operator A. 

The state contributes to the observation y , which also includes a stochastic, 
additive measurement noise v: 



y[k] =C-x[fc]+v[A;]. 



(2) 



In the classic Kalman filter the process and measurement noises are assumed to be 
Normal processes with known co- variances W and V. 1 

Now suppose that we have an estimate x[fc — 1] of the state, and also an estimate 
of the error co- variance P[k — 1] in the estimate, at sample k — 1. The Kalman 
filter uses these estimates, the observation y[k] at sample k, and A, C, W and V 
to form an estimate of the state and its error co-variance at sample k: 



where 



x[fc] := K[k] ■ (y[k\ - y[k\) 

P[k] := (I - K[k] ■ C) • P[k] ■ (I - K[k] ■ Cf 



(3) 
(4) 



y[*] 

K[k] 
P[k] 



C • A • x[k - 1] 

P[k]-C T / (V + C-P>]C T ) 

A-P[fc-1] • A T + W. 



(5) 
(6) 
(7) 



The estimated system state x[/c] (e.g., the generalized coordinate and conjugate mo- 
mentum of the violin mode normal mode of the wire) is thus completely determined 
by the observation y[k], the estimated state at sample k — 1, the wire dynamics, 



^ Even when w and v are not Normal the Kalman filter estimates of the state x[fc] can be shown 
to have the smallest mean-square error of all linear state estimators that depend only on the 
co- variances. 



and the statistical properties of the process and measurement noise. The error in 
the estimate x[/c] falls with k, converging upon a limiting error covariance that is 
fully determined by A, C, W and V; correspondingly, we can choose any initial 
estimate of x and P and the filter will, after several iterations, adjust the state 
estimate and error accordingly. 

From the state estimate at each sample we can, through the measurement equa- 
tion, estimate the contribution of the system to the actual observation. This es- 
timated contribution can be subtractively removed from the actual observation, 
leaving a residual that is as free from the contaminating influence of the process as 
we can make it. 



Ill MODELING THE VIOLIN MODES 

To describe a Kalman filter for our system we need a model for the wire state 
dynamics, the relationship between the wire state and the appearance of the mode 
in the detector data channel, and estimates of the process and measurement noise. 

The wire motion contributes significantly to the data channel only in a narrow 
band about the violin mode resonant frequency; correspondingly, we can focus 
attention this narrow band and model the wire dynamics as a viscously damped 
harmonic oscillator driven by white noise 2 : 

$ = -w&->ri> + N(t). (8) 

We assume that the measurement is of the state variable ip plus white measurement 
noise. 

There are many violin modes, corresponding to the many wires that are used to 
suspend the interferometer mirrors: for each wire there are separate state variables 
and a separate equation describing the dynamics of that mode. 

Using standard lock-in techniques we mix the ip and the data channel with a 
local oscillator whose frequency is near that of the violin modes and band-limit the 
output of the lock-in to the narrow band over which our model is accurate. The 
in-phase and quadrature-phase components of mixed-down ijj become the state x 
used in the Kalman filter, and the in-phase and quadrature-phase components of 
the data channel become the observation y. The matrices A and C, describing 
the state dynamics and the relationship between the state x and observation y are 
derived from the our model in light of the lock-in, band-limiting, and discrete time 
sampling operations. 



2 ) Over the relevant, narrow band near the resonant peak there is no distinction between viscous 
and structural damping. 



IV RESULTS AND DISCUSSION 



A Removing the artifact 



To explore the effectiveness of the Kalman filter in identifying the contribution of 
the violin modes to the detector output we have applied it to data taken in Novem- 
ber 1994 at the LIGO 40M prototype detector. In this instrument the fundamental 
violin mode resonances are all in the (571.6, 605.425) Hz band. 

The upper and lower panels of figure 1 show the power spectra of the interfer- 
ometer output in a 45 Hz band between 565.0 and 610.0 Hz before and after the 
subtractive removal of the Kalman filter estimate of the violin mode contribution 
to the detector output. The filter identifies the contribution of the mode to the de- 
tector output, allowing us to suppress this artifact by 40 dB. The residual bumps 
positioned in the wings of the removed lines, are non-linear artifacts: the violin 
mode amplitudes are so large that they modulate the detector transfer function, 
up-converting other detector noise frequency-modulating the violin mode signal 
itself. 



Caltech 40m data, 19nov94.4/Locked004.nc, 20.43 minutes, f =9868.421 Hz. 
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FIGURE 1. Power spectrum of the 40 meter data showing the violin modes between 571 and 
605 Hz before and after subtraction of the Kalman estimates for all the modes. The remaining 
line feature is the 9th harmonic of the 60 Hz power main. Features like these are dealt with in 
other ways [2-5]. 



B Statistics of the artifact and residual 



Some simple exploratory statistics show the value of identifying and removing 
the known instrumental artifacts from the data stream. 

Figure 2 shows a histogram of the sample amplitude relative to the RMS sam- 
ple amplitude for the data channel before (top) and after (middle) removal of the 
Kalman filter estimate of the violin mode contribution, and (bottom) for the esti- 
mated violin mode contribution itself. In an ideal world the measurement noise and 
the process noise are Normal; correspondingly, each of these distributions should be 
Rayleigh and appear, in these figures, as straight lines. Departures from a straight 
line thus imply non-Gaussian noise statistics. 

Comparing the three panels in figure 2 shows that the violin mode artifact con- 
tributes significantly to the non-Gaussian component of the noise in the detector 
data channel. Since gravitational waves do not excite the violin modes, this excess 
noise component is strictly technical. Generalizing to the full-scale detector, remov- 
ing artifacts like these, together with their associated excess noise, from the data 
stream before analysis thus strengthens our ability to make significant statements 
regarding the detection of gravitational waves in the residual. 



FIGURE 2. Histograms of the square of Kalman estimated amplitude divided by mean am- 
plitude for (top) data in the narrow-band (570-599.5 Hz) having the fundamental violin modes, 
(middle) residual after subtraction of the Kalman estimates of the modes from the data and (bot- 
tom) the Kalman estimated violin modes. Deviations from Gaussianity occur with much lower 
frequency in the residual than in the mode contribution to the observation. 



Kalman estimator and residual statistics, 19nov94.4/Locked004.nc,20.43 minutes , t =9868.421 Hz. 




C Violin modes vs. gravitational waves 



The Kalman filter estimates the violin mode contribution to the gravitational 
wave channel. If that estimate is influenced by the presence of a gravitational wave 



signal then removing the estimated contribution may distort evidence of the wave 
in the output. 

The Kalman filter identifies the violin mode contribution through its dynamics. 
Since the evolution with time of expected gravitational wave signals is different 
than the dynamics of the mode contribution to the detector output, we expect that 
the Kalman estimates of the violin mode contribution will not be influenced by the 
presence of a gravitational wave signal. 

It is useful to consider two different kinds of sources: burst sources, such as 
inspiraling neutron star binary systems, and periodic sources, such as a pulsar at 
a frequency near to but not identical with the violin mode frequency. 

Nearly all the signal to noise ratio signal-to-noise from an inspiraling binary 
is contributed when the signal is in the band from about 70 Hz to 250 Hz [6,7]: 
what happens in the band near the violin mode is is inconsequential. Similarly, 
the filter output is entirely unaffected by what happens in the band where the 
S/N is deposited. We have verified this by forming estimates of the violin mode 
state from the LIGO 40M prototype data, and from the same data set but with an 
added, simulated gravitational wave signal corresponding to a coalescing neutron 
star binary. Even for a very strong signal there is no difference in the predicted 
mode amplitudes when the signal was in the relevant band. 

In separate experiments we have looked at how well the Kalman filter rejects 
nearby monochromatic signals, such as might arise from a pulsar. As long as the 
signal is greater than a linewidth away from the violin mode, the estimated state is 
unaffected by the periodic signal; correspondingly, the periodic signal is unaffected 
the subtractive removal of the estimated contribution of the mode to the output. 

D Monitoring the wire state 

The Kalman filter estimates separately the state of each violin mode; correspond- 
ingly, we can monitor the state of each of these wires separately. It is convenient to 
represent the state in terms of its amplitude and phase (as opposed to generalized 
coordinate and momentum); additionally, it is convenient to express the amplitude 
as an instantaneous measure of the energy in the mode, expressed in terms of tem- 
perature. Excess noise will raise the effective temperature of the mode, while the 
projection of the corresonding mirror motion normal to the optical axis will lower 
the effective temperature. 

V SUMMARY 

Thermal and technical noise that disturbs the violin modes of the substrate 
suspension moves the mirrors, leading to a strong, narrow band artifact in the 
gravitational wave channel. Gravitational waves, on the other hand, also change 
the distance between the mirrors, but without moving or otherwise exciting these 
modes. The Kalman filter described here distinguishes between excitations due 



to gravitational waves and those due to thermal or other excitations of the violin 
modes, allowing us to eliminate a broad class of technical noise sources that might 
otherwise masquerade as gravitational wave bursts. 

A Kalman filter uses the known dynamics of the modes to distinguish between 
the mode "signal" and other contributions to the measured detector output: i.e., it 
detects the violin modes. This distinguishes it from other methods (e.g., multitaper 
methods, linear notch filters [8]) which purport to characterize or remove artifacts, 
but which in fact simply suppress all contributions to the noise within a narrow 
band, not distinguishing violin mode from other contributions. 

The computational cost of identifying and removing the violin modes using the 
Kalman filter described here is negligible: an interpreted Matlab [9] implementa- 
tion on a low-end workstation runs at greater than 20 x the detector's real-time 
sample rate. A compiled implementation, with attention paid to optimization, an 
additional speed-up of 10 or more can be expected. 
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